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We study multifield aspects of Dirac-Born-Infeld (DBI) inflation. More specifically, we 
consider an inflationary phase driven by the radial motion of a D-brane in a conical throat 
and determine how the D-brane fluctuations in the angular directions can be converted into 
curvature perturbations when the tachyonic instability arises at the end of inflation. The 
simultaneous presence of multiple fields and non-standard kinetic terms gives both local and 
equilateral shapes for non-Gaussianities in the bispectrum. We also study the trispectrum, 
pointing out that it acquires a particular momentum dependent component whose amplitude 
is given by /^£/^x- We show that this relation is valid in every multifield DBI model, 
in particular for any brane trajectory, and thus constitutes an interesting observational 
signature of such scenarios. 
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I. INTRODUCTION 

Current measurements of the cosmic microwave background (CMB) anisotropies, such as those 
obtained by the WMAP satellite, already provide us a with a wealth of valuable information about 
the very early universe. Furthermore, with the successful launch of the Planck satellite and the 
increasing precision of large scale structure surveys, one can hope to get yet more precise infor- 
mation in the near future. In this context non-Gaussianity [1] is particularly exciting since it has 
the ability to discriminate between models which are otherwise degenerate at the linear level: for 
instance, a detection of so-called local non-Gaussianity would rule out all single field scenarios of 
inflation in a model independent way [2]. Amongst such single field models, those with standard 
kinetic terms, in which the inflaton slowly rolls down its potential, come with an unobservably low 
level of non-Gaussianity and hence are still consistent with present-day observations j3j. Theoret- 
ically though, their embedding in high-energy physics theories is hampered by the eta-problem, 
namely that Planck-suppressed corrections often lead to potentials that are too steep to support 
slow-roll inflation. In the context of string theory, this problem was demonstrated to be particu- 
larly acute in slow-roll inflationary models based on the dynamics of D-branes moving in higher 
dimensional spaces |H E|. Nonetheless, this set-up precisely motivates an interesting way to bypass 
the eta-problem that has attracted a lot of attention: since the inflaton field in brane inflation is 
governed by a Dirac-Born-Infeld (DBI) action characterized by non-canonical kinetic terms, there 
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exists an upper bound on the inflaton velocity that allows one to achieve an inflationary phase with 
otherwise too steep potentials [6]. Besides alleviating the eta-problem, the non standard kinetic 
terms also enhance the self-interactions of the inflaton, resulting in significant non-Gaussianities 
of equilateral type [7] . Indeed, the simplest single-field DBI models are already under strain from 
observations [8] [9] . 

However, as is central to the discussion of this paper, DBI inflation is naturally a multifield 
scenario, since the position of the brane in each compact extra dimension gives rise to a scalar 
field from the effective four-dimensional point of view \10\ [TTJ (T2J [13]. In multiple field models, 
the scalar perturbations can be decomposed into (instantaneous) adiabatic and entropy modes by 
projecting, respectively, parallel and perpendicular to the background trajectory in field space [14 . . 
If the entropy fields are light enough to be quantum mechanically excited during inflation, they 
develop super-Hubble fluctuations that can be transferred to the adiabatic mode on large scales. 
This effect, as well as the nonlinearities imprinted on the fields at the epoch of horizon crossing, 
was taken into account in [15] where it was shown that the shape of equilateral non-Gaussianities 
is the same as in the single-field case while their amplitude is reduced by the entropy to curvature 
transfer, which therefore eases the confrontation with the data. This paper, as well as subsequent 
ones on multifield DBI inflation [TBI [T7] dS] [TU] EU] EJ [22], focused mainly on equilateral non- 
Gaussianities. More generally, multiple-field inflationary models are known to produce possibly 
large non-Gaussianities of another shape, namely local non-Gaussianities, that arise due to the 
nonlinear classical evolution of perturbations on super horizon scales. The formalism developed in 
these papers also remained very general and no particular model was presented for the entropic 
transfer on large scales. We address both these questions here, using a mechanism outlined by 
Lyth and Riotto [23j to convert entropic into adiabatic perturbations at the end of brane inflation. 

In this scenario, inflation is still driven by a single inflaton scalar field, namely the D3-brane 
solely moves along the standard radial direction of the throat in the context of warped conical 
compactifications. When the mobile D3-brane and an anti D3-brane sitting at the tip of the throat 
come within a string length, an open string mode stretched between them becomes tachyonic, 
triggering their annihilation and the end of inflation [24 . As the brane-antibrane distance is six- 
dimensional, it acquires some dependence upon the fluctuations of the light fields parametrizing 
the angular position of the brane. Hence the value of the inflaton at which the instability signals 
is modulated and the duration of inflation varies from one super-Hubble region to another. In this 
way initially entropic perturbations are converted into the curvature perturbation. The relevance 
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of this mechanism in the slow-roll 'Delicate Universe' scenario jl] has been recently investigated 
|25j . Here however we would like to combine it with the DBI inflationary regime. This was looked 
at by Leblond and Shandera [26 regarding the power spectrum. In the present paper, we extend 
their analysis by taking into account the enhancement of the angular fluctuations by the low speed 
of sound [15 , as well as by investigating the non-Gaussian properties of the curvature perturbation. 
It should be noted that while no explicit model of DBI inflation in a string theory framework has 
been constructed yet that both satisfy observational constraints and are theoretically self-consistent 
|105] - mainly because of the existence of a geometrical limit for the size of the throat in Planckian 
units \27\ [2H] |lU6j - it was so far assumed that the fluctuations of the primary inflaton create 
the curvature perturbation. One can therefore hope to embed consistently the DBI inflationary 
scenario in string theory with other mechanisms to generate the density perturbations, such as the 
one considered in this paper. 

In our model, the curvature perturbation is nonlinearly related to the entropy perturbations, 
therefore the conversion process gives rise to local non-Gaussianities besides the standard equilat- 
eral ones generated at horizon crossing. To the best of our knowledge it is the first time that definite 
predictions are made for an inflationary scenario where both significant local and equilateral non- 
Gaussianities can arise, characterized in the bispectrum respectively by the parameters and 
ftfL' While the logical possibility that a linear combination of different shapes of the bispectrum 
can arise is an acknowledged fact, it should be stressed that what we consider is not merely a 
juxtaposition of an inflaton with non standard kinetic terms generating f^ L and another light 
scalar generating f l £ c L . In multifield DBI inflation, light scalar fields other than the inflaton and 
with derivative interactions naturally contribute to both types of non-Gaussianities. We show that 
this nontrivial combination leaves a distinct imprint on the primordial trispectrum, which acquires 
a particular momentum dependent component whose amplitude is given by the product f^r ■ 
This relation is structural and is valid independently of the details of the inflationary scenario, i.e. 
for any brane trajectory and any process by which entropic perturbations feed the adiabatic ones. 
Hence it constitutes an interesting observational signature of multifield DBI inflation. 

The layout of this paper is the following. In section |Hj we describe our set-up and recall results 
concerning the amplification of quantum fluctuations in multifield DBI inflation. We explain the 
mechanism by which entropic perturbations are converted into the curvature perturbation at the 
end of brane inflation. Using the 5N formalism, we also derive the relevant formulae for quantifying 



this effect. In section III , we calculate the power spectrum of the primordial curvature perturbation. 



We show, in particular, that the entropic transfer is more efficient in the DBI than in the slow- 
roll regime. We also calculate the primordial bispectrum which acquires a linear combination of 
both the local and equilateral shapes of non-Gaussianities. We then combine our results for the 



spectrum and bispectrum to derive constraints on the model. In section IV we turn to the study 
of the primordial trispectrum. We calculate the local trispectrum parameters ttvl and g^L and 
discuss the purely quantum contribution coming from the field trispectra. Then we point out the 
presence of a particular component of the trispectrum whose amplitude is given by the product 
JNL /atl - We finish by plotting the corresponding shape of the trispectrum in different limits, and 
this turns out to have characteristic features. We summarize our main results in the last section. 



II. GENERATING THE CURVATURE PERTURBATION AT THE END OF BRANE 

INFLATION 

A. The set-up and the amplification of quantum fluctuations 

Our setting is that of a flux compactification of type IIB string theory to four dimensions [29 , 
resulting in a warped geometry in which the six-dimensional Calabi-Yau manifold contains one or 
more throats. The ten-dimensional metric inside a throat has the generic form 

ds 2 = h\y K ) 9flu dxV + h-\y K ) G u (y K ) dy 1 dy J , (1) 

where g^ v is the metric of the four-dimensional, non-compact, spacetime and we have factored 
out the so-called warp-factor hdj 1 ) from the metric Gjj in the six compact extra dimensions. In 
the following, we assume that the geometry presents a special radial direction, in agreement with 
known solutions of the supergravity equations [30] , so that the warp factor is a function of a radial 
coordinate p only, decreasing along the throat down to its tip at p = [107 . In this framework, 
we consider the following internal metric 

Gu(y K ) dy'dy 1 = dp 2 + b 2 (p) ggld^d^ , (2) 

where we refer to b(p) as the radius of the throat and ip m (m = 5,6,7,8,9) denote its angular 
coordinates. Since we aim to present a general mechanism, we do not specify a precise form for 
h{p) and b{p) and only require that they approach constant values /iti p and fotip near the tip, which 
is the situation encountered for instance in the Klebanov-Strassler throat [30 . 
We take then a probe D3-brane, of tension 
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where m s is the string mass and g s the string coupling, filling the four-dimensional spacetime, and 
point like in the six extra dimensions. The D3-brane has coordinates y,^ and falls down to the tip 
of the throat where a static D3 is sitting. Following [15J, in terms of the rescaled scalar fields 

¥ = y/nvfa - <t> = VT3 P{b) ,e m = ^n^ b) , (4) 

and rescaled warp factor 

f=(T 3 h±)- 1 , (5) 

the D3-brane low-energy dynamics is captured by the Lagrangian 

P = -fW 1 (Jtetffl + fG I jd^ I d u <pJ) -lj- V^ 1 ) , (6) 

where V(4> 1 ) is the field interaction potential. Note that in general, there are also contributions 
from the presence of various p- forms in the bulk as well as the gauge field confined on the brane. 
In [18], these fields were shown to have no observable effects on scalar cosmological perturbations 
at least to second-order so we have omitted them here for simplicity. 

The explicit calculation of the potential V in ^ is extremely difficult and requires a detailed 
knowledge of the compactification scheme (see e.g [31] and references therein). In general though, 
we know that bulk as well as moduli stabilizing effects break the isometries of the throat, stabiliz- 
ing some of the angular coordinates of the branes. However, there typically remain approximate 
residual isometries of the potential, as shown explicitly in |25j for the most-advanced brane in- 
flation model [5]. For simplicity, we consider only one such isometry direction ip, entering the 
five-dimensional metric glnn through 

g^ m dr = d^ + ... (7) 

and we assume that the four other brane angular degrees of freedom are frozen in their minima of 
their effective potential at the position of the antibrane along these directions. Therefore we take 
the potential V = V{4>) to depend only on the radial position of the brane, which itself moves along 
the radial direction only - eft ^ , 9 m = 0. As opposed to single-field inflation, the perturbations 
along the isometric, i.e. flat direction ip, can be quantum mechanically excited during inflation. In 
that case, the angular D3 — D3 separation 9 = ^JT^/S.'^p varies from one Hubble patch to the others. 
Although this does not modify the dynamics during inflation, this will turn out to be crucial at 
the end of inflation, as we will see below. 
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We now recall the relevant results of [16J regarding the amplification of quantum fluctuations 
in multifield brane inflation, in particular the amplitude of the inflaton and of the angular per- 
turbations at horizon crossing. For that purpose, it is convenient, after going to conformal time 
r = J dt/a(t), where a(t) is the cosmological scale factor, to work in terms of the canonically 
normalized fields given by 

Va = Q<t>> Vs = ~7= Qs ' ( 8 ) 
c s ' V c s 

where Q denotes the perturbations of the field <fr in the flat gauge and 

C S = yJl-fWft (9) 



is the propagation speed of scalar perturbations (see Eq. (10) below), or speed of sound. Clearly, 



from Eq. (9), there is an upper bound on the inflaton velocity \(p\ < -j== . When c 2 ~ 1, one 
can expand the square-root in the Lagrangian Q to quadratic order in the fields. Then the action 
becomes canonical and one recovers the slow-roll regime. However, when the brane almost saturates 
its speed limit — <? s <C 1 — the non-standard structure of the action <JgJ) must be fully taken into 
account: this is the relativistic, or DBI regime we are particularly interested in. 

In Fourier space the equations of motion for v a and v s at linear order take the simple form 

HUGS] 

v" a + (c 2 s k* - ~) v a = , < + (c^fc 2 - Z l + a* A v s = 0, (10) 

where we have introduced the two background-dependent functions z(t) = a^/2e/c s , z s (t) = 
a/y/c^, with e = —jp the inflationary deceleration parameter. The effective entropic mass squared 
H 2 S is given by 

where b(4>) is the radius of the throat evaluated at the brane position. In the following we assume 
that the time evolution of e and c s is very slow with respect to that of the scale factor [108], as 
quantified by the slow-varying parameters 

1 = ^« 1, (12) 
. = (13) 

so that z" / z ~ z's/ z s — 2/r 2 . The amplification of the vacuum fluctuations at horizon crossing 
is possible only for very light degrees of freedom. Although this is automatically verified for the 
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adiabatic perturbation v a because of our assumption z" j z ~ 2/r 2 , this is not necessary true for 
v s \ if nl is larger than H 2 , this amplification is suppressed and there is no production of entropy 



modes. Note that the effective entropic mass squared (11 ) is non zero despite the angular direction 



being exactly isometric. In particular, because the potential and radius of the throat typically 



increase and / typically decreases with cf>, the first two terms in (11 ) are positive. Below we assume 
that \fj, 2 \/H 2 <C 1 so that the entropy modes are effectively amplified. 

Following the standard procedure (see e.g. [33] or [E]), one then selects the positive frequency 



solutions of Eq. (10), which correspond to the Minkowski-like vacuum on very small scales: 

""^vfe^ 1 -^)- <14) 

As a consequence, the power spectra for v a and v s after sound horizon crossing have the same 
amplitude. However, in terms of the initial field perturbations, one finds, using Q, 

(the subscript * indicates that the corresponding quantity is evaluated at sound horizon crossing 
kc s = aH). Therefore, for small c s *, the entropic modes are amplified with respect to the adiabatic 



modes. As we will discuss in subsection II C the standard formulae of the 5N formalism are 
expressed in terms of fields whose perturbations have the canonical amplitude H*/2n at horizon 
crossing. We therefore define the 'canonical' entropy field as 

E = b*c s *6 . (16) 



B. The conversion process 



We now address the question of how initially entropic perturbations can be converted into the 
adiabatic modes. We assume that inflation does not end not by the breakdown of the slow-roll 
conditions but rather persists all along down the throat. Then, when the D3-brane comes within 
a string length of the anti D3-brane, a tachyonic instability arises which ends inflation. Using 
Eqs. ([I]), ([2]) and Q, this happens when 

rj- ((Ap) 2 + 6 t 2 ip (AV) 2 ) = l 2 s , (17) 

n tip 

where A.p and are the radial and angular D3 — D3 separation and l s = m~ l is the string length. 
In terms of the rescaled fields, the tachyon surface, represented in Figs. ([TJ and ([2]), is given by 
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FIG. 1: A simplified picture of the geometry at the tip of the throat, with one angular direction only. 
The radial inflationary trajectory is represented by the blue line. 

\m 

2 + b 2 ip 6 2 = T 3 l 2 h 2 ip . (18) 

The main point is that the end value of the innaton acquires a spatial dependence through the 
fluctuations of the light angular field 6 (see Fig. ([2])). Consequently, the duration of inflation varies 
from one super-Hubble region to another and this can be interpreted as a curvature perturbation, 
as we will quantify in the next subsection. To simplify the notation, we define 

4> C = V^Utip = mskt ' 1P 1/2 (19) 

(2^)3/2 5 y 2 

and the angle —tt/2 < a < tt/2 such that the background value of the innaton field when it reaches 
the tachyon surface <f> e is given by 

4> e = cos(a) 4> c . (20) 

(here and in the following, the subscript e indicates the end of infation). For example, the D3 — D3 
angular separation vanishes when a = 0. 

Let us now make some remarks regarding the validity of our scenario: clearly the D3-brane 
will reach the tachyon surface during its radial fall-down only if the background angular brane 
separation Aip verifies 

A^< l -p^. (21) 

Otip 
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FIG. 2: The tachyon surface, at which inflation ends, in the (f> — 6 plane. The end value of the inflaton is 
shifted from (f> e to <p e + 8<j) e due to the angular fluctuation Qg, hence the duration of inflation varies from 
one supcr-Hubble region to another. 



For the supergravity approximation to be valid, the radius at the tip bf ip must be large in local 
string units ht\ v l s (it is ~ \Jg s M in the KS throat with M>la flux integer). Therefore, condition 



Eq. (21) is a non-trivial requirement, which we assumed to be fulfilled. Note however that if this 
is not satisfied, inflation could end through the angular motion of the brane at the tip |35[ l36l 137] . 
The angular fluctuations would then still produce a spatially-dependent time-delay to the end of 
inflation, and one could expect similar effects to the ones discussed below to arise. Finally, we have 
also assumed for simplicity that the trajectory is completely radial until the end of inflation: more 
generally the trajectory can be nontrivial in the angular directions, for instance if the Coulombic 
attraction between the branes becomes important towards the end of inflation. In that case, 
the entropic perturbations would feed the curvature perturbation through the bending of the 
inflationary trajectory |14] [58] . This mechanism to convert entropic into adiabatic perturbations 
is different from the one used in this paper but both can be present, as in multi-brid inflation 

[391 SQl SD S2] [HQ] 

C. 5N formulae 



In order to quantify the curvature perturbation generated by the entropic fluctuations, it is 
convenient to use the 5N formalism [33] 2H 23 23 2Z] 5 in which a key role is played by the local 
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integrated expansion, or local number of e-folds, between some initial and final hypersurfaces 



AT(x) 



H(t,x)dt. 



(22) 



In this formalism, the curvature perturbation on uniform energy density hypersurfaces, which we 
denote £, is identified as the perturbation in the local number of e-folds of expansion from an 
initially flat hypersurface to a final uniform energy density hypersurface 



( = 5N = JV(x) - N , 



(23) 



where N is the number of e-folds in the homogeneous background spacetime. In the long- wavelength 
limit, according to the separate universe picture [48 , the integrated expansion can be calculated 
from solutions to the unperturbed Friedmann equation, with initial conditions specified by the 
perturbed scalar fields ( = N(tp A ) — N, where (p A = cp A + Q A is the sum of the homogeneous values 
plus fluctuations of the scalar fields on the initial spatially-flat slice, which we take to be soon after 
horizon crossing during inflation. Taylor-expanding this relation in terms of the field fluctuations 
leads to the formal expression 



C = N A Q A + l -N AB Q A Q B + \n abc Q a Q b Q c + ... 
2 6 



(24) 



Once N(ip A ) is known, one can work out the coefficients N A , Nab, N A bc in (24) and determine 
the curvature perturbation. 

In our model, since the background dynamics is solely determined by the radial scalar field 4>, 
the local number of e-folds, evaluated right after the end of inflation, simply reads 



AT(0*(x),H*(x)) 



(25) 



where the canonical entropy field S was defined in (16) and, using Eqs. (18) and (19), the end 
value of the inflaton reads 



i— 'si 



with 



= 



9tip 



(26) 



(27) 



being the ratio between the radius at the tip of the throat and at sound horizon crossing. Here 
we assumed that no significant expansion is generated after the field reaches the tachyon surface, 
that is, the sudden end approximation [49 . If not negligible, the extra curvature perturbation 
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generated can be taken into account similarly to [39], though even at this stage, £ may not have 
settled down to its final value. For example, the further evolution of the tachyon can give other 
contributions to £ JSOJ EI] . This interesting aspect is outside the scope of this paper but it should 
be borne in mind that it is present in general. 



From Eq. ( 25 ) , one obtains 



where 

c* = 

and 



4>*+Q<t>* /ff 



H 



C = C* + Ce 



1 d (E 

2d0 



q1* 



i d 2 



3! dcp 2 V 4> 



H 



(28) 



QL + — (29) 



e(H + Q H ») f H 
0e(H) 



H 



+ 



1 d ( H 

2d0 



,2 1 d 2 



3! 



+ ... 



(30) 



are the contributions to the curvature perturbation from respectively, the epoch of horizon crossing 
and the end of inflation. Loosely speaking, £ e is associated to the time delay generated by the 
fluctuation 8(j) e in Fig. where we use the notation 



(31) 



and the derivatives of </> e (H*) (26), given in Appendix 1, are evaluated on the background. Note 
that in order to trust the perturbative expansion, the angular fluctuations must be small compared 
to the background separation, namely j3H^/c s ^ <C 4>c- 



Adding (29) and (30) gives an expression of C of the form (24), where A, B = a, s with Q a = 



b t and Q s = Qs,, which are normalized to share the canonical amplitude H*/2ir. To leading 



order in the slow-varying parameters r] and s (12)-(13) and their time derivatives, only the first 



terms remain in the expansions (29) and (|30|) so that the non-zero coefficients in (|24|) are then 

1 



N ss 



1 



V2eCsM pl 

A 2 ) 



1 



J?) 



(32) 
(33) 



where we have used H/(j> = —(2ec s )~ l l 2 /M p \ [16]. For completeness, we include in Appendix 1 
the full expansion of Q, not restricting to leading order in the slow-varying parameters. Notice 
also that they cannot be neglected in the computation of the scalar spectral index n s and running 
non-Gaussian parameter rtNG as they give then the leading order result (see below). 
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III. POWER SPECTRUM AND PRIMORDIAL NON-GAUSSIANITIES FROM THE 

BISPECTRUM 



According to Eq. (24), the statistical properties of the curvature perturbation £ are inherited 
from those of the field fluctuations Q A at horizon crossing. For clarity, we first recall the results, 
determined in |16j . for the field two-point and three-point functions in two-field DBI inflation, 
before calculating the power spectrum and bispectrum of C in the following subsections. The 



investigation of the trispectrum will be the subject of section IV We use the notations of [52] . 



A. Statistical properties of the field perturbations at horizon crossing 

In Fourier space the power spectrum of the scalar field perturbations is defined by 

(QkQk') = C AB (k)(2n) 3 5 3 (k + k') . (34) 

To leading order in the field perturbations 

C AB {k) = ^ 3 5 AB , (35) 

where 5 AB is the Kronecker delta-function. Notice that the cross-correlation between adiabatic and 
entropy modes is zero for the straight line background trajectory considered here as the coupling 



between them exactly vanishes in that case (see Eq. (10) as well as [16] for details). 



The bispectrum of the field perturbations is defined by 

(Q& Qi 2 Q% 3 ) = B ABC (h, k 2 , A :3 )(2vr) 3 5 3 (k 1 + k 2 + k 3 ) . (36) 

In slow-roll inflation, where the self interactions of the fields are suppressed by the flatness of the 
potential, the bispectrum of the fields is small, both for single |53j and multifield inflation |54J. On 
the contrary, self-interactions are enhanced in models with non-standard kinetic terms [55], and in 
DBI inflation in the low sound speed limit in particular, with the result [16] 



B ABC (h,k 2 ,k 3 ) = b ABC (h,k 2 ,k 3 ). (37) 

The fully adiabatic momentum dependent factor b ABC (k\, k 2 , k 3 ) is given by 

b^(k u k 2 , k 3 ) = - 1^ [6fc?*|fc§ - fcl(k! • k 2 )(2k t k 2 - k 3 K + 2K 2 ) + perm.] (38) 

where K = k\ + k 2 + k 3 and the 'perm.' indicate two other terms with the same structure as the 
last term but permutations of indices 1, 2 and 3. Note also that it depends only on the norm of 
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the three wave- vectors as for instance, ki • k2 = \{k 3 — k\ — due to momentum conservation. 
This is the standard result from single field DBI inflation [56 . In the relativistic limit, there exists 
only one other non-zero three-point correlation function at leading order, namely 

b™(k u k 2 , h) = n l 3R3 [2kjk 2 2 kl + kj(k 2 ■ k 3 )(2k 2 k 3 - k x K + 2K 2 ) 

- ki{k\-k 2 ){2k l k 2 -k 3 K + 2K 2 )-kl{k l -k 3 ){2k 1 k 3 -k 2 K + 2K 2 )]. (39) 



B. Power spectrum, scalar spectral index and tensor to scalar ratio 



The power spectrum of the curvature perturbation is defined as 



<CkCk'> =Pc(fc)(27r) 3 5 3 (k + k / ). 



(40) 



The corresponding variance per logarithmic interval in /c-space is given, to leading order in the 
field perturbations, by 



V c {k) 



k 3 



2^'= IS M + 



(41) 



where we have used Eqs. (24) and (35). As in |16| l5?] . it is convenient to introduce the transfer 
function T CTS , such that N s = T as N a . The curvature power-spectrum then takes the form 

"'=pb|( 1+J *-)' (42) 

where T 2 S quantifies the contribution of the entropy modes to the final curvature perturbation. 



This vanishes in single-field DBI inflation while in our case, from Eq. (32), 

T 2 S = — - — tan 2 (a)(3 2 . 



(43) 



Let us comment on the ranges of the various parameters entering Eq. (43). First, since the radius 



of the throat decreases from the UV to the IR end, f3 = btip/K is bounded by one and this bound 
is saturated when the last 60 efolds of inflation take place at the tip of the throat, where b(<p) 
becomes a constant [58]. Second, the entropic transfer depends on the angular D3 — D3 separation 



through the angle a (20). When tan(a) = 0, the angular fluctuations give no time-delay to the 



end of inflation at linear order, as is clear from Fig. and hence the transfer function vanishes. 
In that case, the first effect appears through higher order loop corrections (see Appendix 2). In 
the following, we have in mind that tan(a) = O(l) although we keep formulae general. Finally, 



the result (43) indicates that in slow-roll inflation, where the speed of sound is one, the entropic 



contribution to the curvature power spectrum can be significant, compared to the inflaton one, 
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only if the deceleration parameter e is smaller at the end of inflation than at horizon crossing. 
In DBI inflation, however, the transfer function is amplified by the inverse of the product of the 
sound speed at horizon crossing and at the end of inflation, hence it is more efficient. Note that 



the enhancement of entropic perturbations by the inverse of the sound speed ( 15 ) was crucial in 
deriving this result. 

It is straightforward to calculate the scalar spectral index and tensor to scalar ratio from the 
power spectrum. We find 

and pl| 

r = 16^3^ - 16 ^^^ec se cl , Tl s » 1 (45) 

where the last equality holds in the limit of a large entropic transfer. Hence it is clear that, when the 
curvature perturbation is of entropic origin, the links between the observables and the microscopic 
parameters of the model are completely different from the single field case. 

C. Primordial bispectrum 

1. General definitions 

The bispectrum of the curvature perturbation is defined as 

<Ck! Ck 2 Ck 3 ) = B c {k x ,k 2 , fc 3 )(2vr) 3 ( 5 3 (k 1 + k 2 + k 3 ) , (46) 



where, from Eqs. (24), (34) and (36) and to leading order [59] 
B c (h, k 2 , k 3 ) = NaNbNcB abc (ki, k 2 , k 3 ) 

+ N A N BC N D [C AC (h)C BD (k 2 ) + C AC (k 2 )C BD (k 3 ) + C AC (k 3 )C BD (h)] . (47) 

Observational quantities are usually expressed in terms of the dimensionless non-linearity parameter 
Inl~ generally momentum-dependent - defined by 

B c (h,k 2 , k 3 ) = -f NL [P c (h)P c (k 2 ) + perm.] . (48) 



Hence there are two contributions to fNL'- the first, coming from the first term in Eq. (47), is 



(3) 

related to the three-point functions of the fields at horizon crossing, and we will denote it as f NL - 

,(3) = b NANBNcB^jk^k^) 

Inl 6 (P c (fci)P c (A;2) + perm.) ' 1 ' 
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The second, coming from the second group of terms in Eq. (47), comes from the leading order 



nonlinear relation between the curvature perturbation and the field perturbations, and we will 



denote it as f§ c L (27]: 



floe 
J NL 



5 N A N B N 



AB 



6 (N C N C ) 

(3) 



2 ' 



(50) 



the total fxL being the sum of the two /nl = /jyx + -/jvl- 



Equilateral and local bispectra 



(3) i i 

Let us first discuss fj^ L - if inflation is of slow-roll type when the observables modes cross the 
horizon - <? sif ~ 1 - then fSl is negligibly small. We therefore concentrate on the relativistic regime 
4, < 1, in which case §S7§, (@, (§§ and (0 give p2] 



AS) 

J NL 



(51) 



Here we used the symmetry property b ass {ki,k2,k 3 ) + b sas (ki,k2,k 3 ) + b ssa (ki, k 2 , k%) = 

(3) 

b crcrcr (ki,k2,ks), which implies that the shape dependence of f NL is the same as in single-field 
DBI, as can be understood in a geometrical way [21]. In the equilateral limit k\ = k 2 = k%, 

35 1 1 



f eq 



NL 



108 cL 1 + T2 ' 



(52) 



where T^ s is given in (43) in our model. Hence the entropic to curvature transfer in multifield 



DBI inflation diminishes the amount of equilateral non-Gaussianities with respect to the single 
field case. This comes from the fact that the transfer not only enhances the bispectrum of £ but 
it also enhances its power spectrum by the same amount. Since f^ L is roughly the ratio of the 
three-point function with respect to the square of the power spectrum, ff? L is effectively reduced. 



We now turn to the local shape of the bispectrum: from the definition (50), we obtain (111 



floe 

Jnl 



N*N SS 



(53) 



6 (AT2 + AT2)2 • 

When there is a large entropic transfer - N% 3> N% - (the case where the entropic transfer is small 



floe 
JNL 



is discussed in Appendix 2) Eq. (53) reduces to the single field entropic result 

6 iVf 

where we have used Eqs. (|32l)-(33) in the last equality. Up to the factor J~c^ which diminishes the 



pi 



T » 1 

o~s 



(54) 



effect, this expression is identical to the slow-roll result [23] in the same limit. We have however a 
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concrete model at hand, for which one obtains 



f£i = \ ■ 2 * , J ^^* Tl»l. (55, 
o sm \a) cos(a) V 2 



This shows that f^ c L is always positive, in agreement with the discussion in [60 , and can be sig- 
nificant, given that the inflaton field (p is highly sub-Plan ckian in brane inflation [27] and therefore 

it can even saturate the existing observational 



that 4p! > 1. As we will see in subsection 



HID 



bound and put constraints on the model. 

Let us stress that the bispectrum in our scenario displays a combination of two different shapes 
|112] . First, the classical nonlinear relation between the curvature perturbation and the light 
entropic scalar field gives rise to local non-Gaussianities, that peak for squeezed triangles (k% <C 
/c2 ~ and this is quantified by the parameter f l £ c L - This type of non-Gaussianities constantly 
arises when the curvature perturbation is generated at the end of inflation through light fields other 
than the inflaton (Ml EU [62j E31 IMl E5J ESI E3 ESI ES]- Second, derivative interactions produce 
quantum correlations at the epoch of horizon crossing between modes of comparable wavelengths. 
The associated non-Gaussian signal peaks for equilateral triangles in momentum space {k\ ~ k<i ~ 
k%) and this is quantified by the parameter fjj L . Note that as the local and equilateral signals have 
fairly orthogonal distributions in momentum space |7(J| , observational bounds on each of them can 
be used when they are both present, each one being almost blind to the other |113j . 

3. Running non-Gaussianities 

Besides its amplitude and shape, the scale dependence of the primordial bispectrum is another 
probe of the early universe physics, and recently it has been shown that combining CMB and large 
scale structure observations give interesting constraints on the running of non-Gaussianities |71|I72]. 



In our scenario, while f£ c L (55) is scale-independent, it is clear that fJ L (52) is scale-dependent. 



This can be quantified by the running non-Gaussian parameter, defined as 

mG = dink = - 2s *-rr^ r~ g *~ 2 ^H ' (56) 



where we have used Eqs. (43) and (52) in the second equality. 

As the speed of sound generally decreases with time in models in which the brane goes from 
the UV to the IR end of the throat, s = c s /Hc s < so that ti^g is positive if the entropy modes 
do not feed the curvature perturbation. On the contrary, using the relation rj = 2e — s — -jjj valid 
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in the DBI regime (see [18] for details), one obtains, for a large transfer, 

n NG = -2e*+ 2-^- + -/*-, T 2 S » 1 . (57) 

Hence, if observable modes cross the horizon at the tip of the throat where b and / become constant, 
the running non-Gaussian index rather becomes negative. 

D. Implication of the results 

In the recent paper [73], it was noticed that in some models, brane inflation ends by tachyonic 
instability in the relativistic DBI regime, even if inflation is of slow-roll type when the observable 
modes cross the horizon. In this subsection, we concentrate on this limit, namely c 2 e <C 1, and 
additionally assume that the curvature perturbation is mostly of entropic origin T 2 S 3> 1. Note 
that the deceleration parameter e has been used until now only as a small parameter, quantifying 
how much the inflationary expansion is close to de-Sitter. However, in the DBI regime, it can be 
related to the brane tension and the warp factor. Indeed, since 

H _ <j) 2 

E = ~lp- 2c s i?2 Mp V (58) 
when c 2 = 1 — f4> 2 <C 1, this gives, with Eq. (EL 

6Cs = 2i^ • (59) 



Therefore, from (43) and (59), our hypotheses imply the condition 



2 fT 2 

-^^tan 2 (a)/? 2 »l (60) 

I 3 n tip c s* 

on the parameters of the model. With these assumptions, we now combine the results of the 
previous subsections and discuss the constraints imposed by non-Gaussianities from the bispectrum. 



* NL 

obtain 



From given in (55), together with the definitions of 4> c in Eq. (19) and T3 in Eq. &3«, we 



6 sm z (a) cos(a) H e 

To avoid stringy corrections, one requires that at least htl P ms > 1 [74j . which therefore acts in the 
direction of making observable local non-Gaussianities |114j . This can be made more stringent by 



noting that the power spectrum (42) can be reexpressed, in the limit of a large entropic transfer, 
as 
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Then, Eqs. (61) and (62) together lead to 



fNL ~ 6 2 S i n 3/2(2a) WJ Wj ' ( J 

For instance, taking g s = 0.1 and a = ir/4 - which minimizes the result - together with the 
observed normalization of the power spectrum = 2.41 .10 -9 [3] gives 

^ 3oo (£f (t) 1/2 - < 64 > 



Let us recall that the only hypotheses that we used to derive (64) are that the curvature pertur- 
bation is of entropic origin and that the end of infation takes place in the relativistic regime. It is 
in particular valid for any value of c s * between zero and one. Furthermore, since the Hubble scale 
decreases during inflation, namely H*/H e > 1, then we find the lower bound 

f^° L >300P l/2 . (65) 

Therefore, the WMAP5 observational constraint —9 < < 111 (95% CL) [3 j implies the upper 
bound (3 = 6tip/&* ^ 0.1 in order to avoid too large non-Gaussianities of local type. 



We now turn to equilateral non-Gaussianities, which are diluted by the entropic transfer (52) 

^ * ~°- 3 cik = -°- 3 taJ(a) 5^ ' (66) 



With Eq. (59), this gives 



^-°- 3 ^(t) 2 M^) 4 ' (67) 

One needs a precise model to actually determine the amplitude of f^fr . Let us simply comment 



that if one tunes j3 = b ^ R to a small value to low down (|64|) in the observational range, this 



enhances fxj L - However, in that case, observable modes cross the horizon far from the tip of the 



4 



throat, and one expects a huge hierarchy y^ 2 - J <C 1, which tends to put f^ L within the current 
observational bounds -151 < f^ L < 253 (95% CL) 0. 



IV. PRIMORDIAL NON-GAUSSIANITIES FROM THE TRISPECTRUM 

As the next generation of experiments will be able to probe refined details of the statistics of 
density fluctuations |75[ 176] 177] , the study of the primordial trispectrum is becoming increasingly 

(3) 

important. In DBI inflation, another motivation comes from the fact that f NL acquires the same 
momentum dependence in single- and multiple-field models. Therefore, we cannot observationally 
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differentiate between them with the bispectrum alone. As the degeneracy between models tends 
to be broken as we go to higher-order correlation functions, the investigation of the trispectrum in 
multifield DBI inflation is thus very natural. 

A. General definitions 

The primordial trispectrum is defined as the connected part of the four-point correlation function 
of the curvature perturbation in Fourier space 

(C kl Ck 2 Ck 3 Ck 4 )c = T c (ki, k 2) k 3 , k 4 )(2vr) 3 5 3 (k 1 + k 2 + k 3 + k 4 ) . (68) 

We need twelve real numbers to specify a set of four three-dimensional momenta. However, momen- 
tum conservation, ki + k2 + k 3 + k 4 = 0, eliminates three of them and invariance under rotational 
symmetry eliminate three others. We are thus left with six independent parameters that specify 
inequivalent configurations of the tetrahedron formed by the four k vectors: we will use the set 
{k\, k 2 , &3> k/L, k\2, k 23 } where k\ 2 = |ki + k 2 | and ^23 = |k 2 + ki3 1 . The others kij = |kj + kj| can 
then be reexpressed in terms of them as follows: 



&13 = k 2 4 = J k\ + k\ + k\ + k\ - k\ 2 - k 



2 

23 ' 



ku = k 2 3 , k u = k 12 ■ (69) 

Note that there are geometrical limitations on the parameter space, for instance in the form of 
triangle inequalities (see [78 for details). 

From the general 5N expansion ( |24[ ) , the primordial trispectrum can be evaluated and one finds 

m 

r c (k 1 ,k2,k 3 ,k 4 ) = NaNbNcNdT^v^M^M) 

+ N ab N c N d Ne [C AC (k 1 )B BDE (k 12 ,h,h) + 11 perms] 

+ N ab N CD NeN f [C BD (k 13 )C AE (k 3 )C CF (k 4 ) + 11 perms] 

+ N ABC N D N E N F [C AD (k 2 )C BE (k 3 )C CF (k 4 ) + 3 perms] (70) 

(where we have omitted special configurations of the wavevectors, for instance when any k vector 
or the sum of any two k vectors is zero). Here the T ABCD are the connected four-point correlation 
functions of the field perturbations at the epoch of horizon crossing, defined by 

(Qi Ql QgQ&c = T^ OT (k 1 ,k 2 ,k 3 ,k 4 )(2vr) 3 5 3 (^k i ) . 



21 



Introducing 



and 



_N AB N AC N B N C 
™ L = {N D N°f (71) 



_ 25 N ABC N A N B N C 
9NL = 54 (N D NC)3 (72) 



simplifies the expression of ( 70 ) to 



^(kj.k^ks.kd) = iV A AT B iV c iV D T^ CD (k^^k^ 

+ N AB N C N D N E [C AC {k x )B BDE {k 12 ,k^k A ) + 11 perms] 

+ tjv L [P f (fci3)P c (A3)P f (fc 4 ) + 11 perms] 
54 

+ ^9nl [P f (A&)P c (A&)P f (A*) + 3 perms] . (73) 



In order to compare the various momentum dependences entering (73), we use the form factor 
Tdefined by [78] 

r c (k 1 ,k 2 ,k3,k 4 ) = (2^) 6 P ? 3 (jlyt) T(h,k 2 ,k 3 ,k 4 ,k 12 ,k u ), (74) 

whilst the amplitude of the trispectrum signal is quantified by the estimator t^i for each shape 
component 

1 RT 

— rT(fci, k 2 , &3, fe4, fel2, few) component ► *JVL • (75) 

/C limit 

Here the regular tetrahedron (RT) limit means ki = k 2 = k% = k 4 = k\ 2 = k\& = k \22\ [75] . 

In the DBI regime of brane inflation, the quantum non-Gaussianities of the fields are large - 
as set by their bispectrum B ABC and connected part of their trispectrum T ABCD - and the first 



two lines in Eq. (73) a priori gives sizeable contributions to the primordial trispectrum. Moreover, 



the last two lines in Eq. (73) could be significant as we will see later. Their respective amplitude 



are determined by the two parameters t^l and gNL, which are similar to for the bispectrum 
in that they describe the nonlinearities generated classically outside the horizon. In that sense we 
will occasionally refer to them as the local trispectrum parameters. They contribute to the form 
factor as 

36 

T D ^TNLTiocl + gNL,Tloc2 , (76) 
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where the two local shapes are 



'loci 



9 (k\k\ „ \ 
50^ + llpermS J 



n 



ot:2 



50 
27 

Too 



(77) 
(78) 



i=l 



and the size of the trispectrum for each is 



^NL = 1-18 T~NL , ^NL = 1-08 gJVL • 



(79) 



We study the local and intrinsically quantum contributions to the trispectrum in the next 



subsection while the discussion of the second term in the right hand side of Eq. ( 73 ) is deferred to 
subsection IIV CI 



B. Local and intrinsically quantum trispectra in our model 



We begin by considering the standard local parameters tnl and gNL- From their definitions 



Eqs. (71) and (72), they are given by 



TNL 



N 2 s Nl 



(JV2 + JV?) 8 ' 9NL 54(^2 + ^2)3- 
Notice that using Eqs. ( 104 )-( T06| ) and (32)-(33), the entropic derivatives satisfy 



25 N*N SSS 



3sin 2 (a)iV s 2 s . 



(80) 



(81) 



We therefore obtain 



gNL = sin {a)T NL 
lo 



(82) 



As tnl and gNL are in principle observationally distinguishable [7S], this gives us the possibility 
to deduce the angle a from the observations. Given that g^L < ff^TVL, it should be stressed 
also that one can not obtain a large g^L without having a large tnl, contrary to the cases of the 
curvaton [80], ekpyrotic models [81] or non slow-roll multifield inflation [82 . Specializing to the 



limit in which the curvature perturbation is mostly of entropic origin, one obtains from Eq. (80) 
the single-field entropic result 



TNL 



N 2 

1 * ss 

Nj 



•loc 
NL I > 



tL > l 



(83) 



where f £ L is given in this limit by (55) (and the relation (82) still holds) 
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The intrinsic trispectra of the fields also leave their imprint on the primordial trispectrum, 
through the linear relation between the curvature perturbation and the field perturbations at hori- 



zon crossing (see the first term in Eq. (73)). The calculation of the field trispectra in multifield 
models with non-standard kinetic terms, and in multifield DBI inflation in particular, begun only 
very recently |19[ I2T1 [22] . For example, the authors of |21l [22] considered the trispectra from the 
intrinsic fourth-order contact interaction. However, it was pointed out that there are other impor- 
tant contributions coming from interactions at a distance |19| [78] 183) 184"] . As the full calculation 
of the field bispectra in multifield DBI inflation has not yet been completed, below we simply com- 
ment on the partial results coming from the contact interaction. The important result is that the 
degeneracy between single and multiple-field models is broken at this level, i.e. with our notations 
one obtains 

T f (ki,k 2 ,k3,k4) Nt{T aaaa {\iiMMM) +T2 s (T CTCTSS (k 1 ,k 2 ,k 3 ,k 4 ) + 5 perms) 

+ T^T^!, k 2 ,k 3) k 4 )) 

= N A a {l + Tl) (t™™ (k x , k 2 , k 3 , k 4 ) + TlT ssss (kt, k 2 , k 3 , (84) 



where the last line follows from a particular relation between the mixed term T aass and the adia- 
batic and entropic trispectra, which have the same order of mag nitude ~ N%H$/c% but different 
momentum dependence. Therefore the trispectrum can a priori discriminate between single- and 
multifield DBI inflation. In particular, one can hope that looking at the trispectrum signal in 
different limits of the tetrahedron's parameter space might enable a measurement of T as . As for 



the overall amplitude of these terms though, with = (l^f) 2 + T^ S ) (41 ), this schematically 
gives a contribution to t]y l of the form 

tNL ~ 4^*2(1 + 3* )(1 + aTl) ~ + Tl) 1±^*L ( 8 5) 

where a is a O(l) numerical coefficient. Similarly to the case of the bispectrum, notice that the field 
trispectra having the same amplitude implies that the transfer function T as solely determines the 
multiple- field modification to the single-field result t^L ~ -j- [85 . Therefore, even if one does not 
succeed extracting T^ s from measurements of the trispectrum alone - for example because T% s S> 1 



and the entropic signal (the last term in ( 84 ) ) dominates over the adiabatic one - one can in principle 
combine measurements of the trispectrum and of the bispectrum (52) f^ L m — 0-3 e2 rp+rs ) *° 



determine both c s * and T as . This illustrates how the study of higher-order correlation functions 
can be used to extract information one can not obtain with the power spectrum alone. 
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C. An observational signature of multifield DBI inflation 



Going back to the general expression Eq. (73), one sees that the intrinsic bispectrum of the 



fields also contribute to the primordial trispectrum of the curvature perturbation, i.e. 

T^kzMM) D N AB N C N D N E [C AO {fa)B BDE (k 12 ,fa,fa) + H Perms] . (86) 
Remarkably, this term is solely determined by quantities that already appeared in lower-order 



correlation functions. In particular, while Nab determine the amplitude of (50), the intrinsic 



bispectra of the fields B ABC set the magnitude of f^ L (|49|). Therefore, this contribution can be 



non negligible only if significant classical nonlinearities and quantum non-Gaussianities are present 
in the same model, which explains why it has been neglected before and why we refer to it as 
j,ioc,eq - n k e following. Note that since the shape of the field three-point functions does not 
take a universal form, one cannot extract a number that characterizes the amplitude of y^ oc ' eq 



independently of a model, contrary to tjvl (71), which also depends only on lower-order terms. 



T ■ rp loc, eq j 

In our scenario, 1^ reads 

+N ss N u N 2 s C ss {k l ) (B ssa (k 12 , fa, fa) + B s ™(k 12 , fa, fa)) + 11 perms] . (87) 

Notice that the first two terms, being proportional to N ua , are negligible to leading order in the 
slow-varying parameters and we are left with the second line alone at this order. We have included 
the first line only to stress that the trispectrum breaks the degeneracy between the two shapes 



jyaaa anc j foass en t erm g th e bispectrum, defined respectively in Eqs. (38) and (39). Recall indeed 
that we used the symmetry property 

^(fci, k 2 , fa) + b SCTS (fa, k 2 , fa) + b ssa (ki, k 2 , fa) = b aa(7 (ki,k 2 , fa) (88) 

(3) 

to show that the multifield effects do not modify the shape of f NL compared to the single-field 
case. Here, one sees that because N acr ^ N ss in general, one can not use this identity and the 
primordial trispectrum truly depends on the shape b ass , not on its symmetrized version b aaa . 



Let us now concentrate on the leading-order term, namely the second line in Eq. (87), and write 
its contribution to the form factor as 

(89) 

where 

81 

Tioc, eg = — {fafafaf [b ss °(k 12 , fa, fa) + b sas {k 12 , fa, fa)] + 11 perms . (90) 
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Here, snl is a dimensionless non-linearity parameter that set the amplitude of T\ oc eq , similarly to 



ttvl and g^L for T\ oc \ and T/ OC 2 respectively (see Eq. (76)), and the numerical factor in the definition 
of TJ oc e g is for convenience only. The overal amplitude of this contribution to the trispectrum is 



t l ° c £ eq = 0.09 s NL . 



By noting that (37) can be rewritten in the form 

B ABC (h,k 2; k 3 ) -- 



from (87), one formally obtains 



4cL 



175 AT™ 1 



SNL 



648 ^^(1 + 1^ 



which is exactly the product of n?? (53) and ff? r (52) 



_ floe feq 
SNL — JNL JnL ■ 



(91) 



(92) 



(93) 



(94) 



It is worth emphasizing that, in order to derive the result (94), we solely used the 5N formalism 
together with the general results for the field bispectra jT6j. Hence it is clear that it is valid for 
every DBI model, not restricting to a radial trajectory nor to a specific scenario for the entropy to 
curvature transfer. In particular, we considered the two-field case for simplicity of presentation but 
the proof generalizes straightforwardly to a higher number of light fields. The consistency relation 



(94) is thus structural to multifield DBI inflation. It represents a distinct observational imprint on 
the trispectrum of the presence of light scalar fields, other than the inflaton, and with non standard 
kinetic terms, hence the appearance of both the local and equilateral non-linearity parameters. 



D. Shapes of trispectra 



To test the consistency relation (94), one must observationally disentangle the contribution 
proportional to snl from other components of the trispectrum. The study of its six-dimensional 
parameter space is challenging and in the following, we follow the discussion in [75] and simply 
plot the form factor 7} OCj eq in various limits to reduce the number of variables, and compare it to 
Ti oc i and T/ OC 2. The shape functions are left white when the momenta do not form a tetrahedron. 
We consider the following cases (see |78j for details): 

1. Equilateral limit: k\ = k 2 = k% = k^. In Fig. [3] we plot I] OCj eq , T\ oc \ and T] OC 2 as functions 
of k\2/k\ and fcu/Zci. One observes that Ti oc eq and 7} ocl blow up at all boundaries, cor- 
responding to k\2/ki, kn/ki and ki^/k\ — > 0. However, the signal in T/ OCi eq is much more 
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important, as is clear from the extremely large cutoff for the z-axis in this figure. On the 
contrary, TJ OC 2 is constant in this region of parameter space. 

2. Specialized planar limit: we take k\ = ks = k\4, and the tetrahedron to be a planar quad- 
rangle with [75] 



k%2 



k l + ~^ (k 2 h + ^(4kl-ki)(4kf-kl 



1/2 



(95) 



We plot the shape functions as functions of k%jk\ and k^jk\ in Fig. |4j At the &2 — ► fei limit, 
we have fci3 — > 0, so that I) OC) eq and 7] oc i blow up. Again, the signal in T} OC) eq is much 
more important and we had to use an extremely large cutoff for the z-axis for the sake of 
presentation. The most distinctive difference comes from the fc2 — > and k^ — ► boundaries, 
where T[ oc \ and 7} OC 2 take finite (non-zero) values whereas 7] OCi eq blows up. Moreover, its 
sign alternates non-trivially over the parameter space. Interestingly, this feature was not 
present in any other shape studied in [75] . 

3. Near the double-squeezed limit: we consider the case where fc 3 = k^ = kyi and the tetrahe- 
dron is a planar quadrangle with [75] 



I VI (_u2 , p , t,2\ _ V2 jl2 , ju2 1,2 , u2 u2 , j.2 ju2 _ l2 p _ £.4 , pp 
„ ft l I "'12 ^ ^3 ^ 'MJ fi 'sl"'s2 ^ ^12^14 ^ ft 12"'4 ^ "'14"'4 /i 14"'3 ft 4 ^ ' i 3'M 

h = 1 JE. • (96) 

where k s \ and fe S 2 are defined as 

k 2 sl = 2a/ {MM + ki • k i )(kik 4: - ki • k 4 ) , 

fc s 2 2 = 2A/(fc 3 /c 4 + k 3 • k 4 )(k 3 k 4 - k 3 • k 4 ) . (97) 

We plot Ti OCi e q , Ti oc i and T\ OC 2 as functions of k 4 jk\ and k\ 4 jk\ in Fig. [5} Similarly to 
the specialized planar limit, notice that the sign of 7] OC) eq varies non trivially over the 
domain. Several differences between the shapes are visible. 1) In the squeezed limit, at 
(ki/ki = 1, k\/±jk\ = 1) where k2 — ► 0, and in the double-squeezed limit, k% = k 4 — > 0, 7} OCj eg 
blows up while the local shapes are finite. 2) In the folded limit {k^jk\ = \,k\t±jk\ = 0), 
both 7] OCi e „ and T\ oc \ blow up and T\ OC 2 takes a finite value. 3) In the other folded limit, 
{k/±jk\ = l,ki 4 /ki = 2), all shapes remain finite. Close to it, note that the bump in 7] oc i 
actually remains finite while T] OCi eq truly blows up. 

To summarize, although we did not make an exhaustive study of the shape factor 7} OCi eq , we have 
seen that it presents characteristic features. In particular, in the last two envisaged situations where 
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FIG. 3: In this group of figures, we consider the equilateral limit ^1 = ^2 = ^3 = k^, and plot Ti oc ^ eq , 
Tied and Ti OC 2, respectively, as functions of k\ijk\ and ku/kx. Note that T; OCj eq and J] ci blow up when 
ki2 <C k\ and k\/± <C fci, as well as in the other boundary, corresponding to k\^ <gC k\. 
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FIG. 4: In this group of figures, we consider the specialized planar limit with k\ = k^ = kn, and plot 
Tioc, eq, Tiod and T/ OC 2, respectively, as functions of kilk\ and ki/ki. Along the diagonal k2 — > Ti OCt eq 
and Tiod blow up because in this limit, fci3 — + 0. At the &2 — > and — > boundaries, TJ OC] eq blow up 
while the local shapes remain finite. Notice that the sign of T\ oc eq varies non trivially over the parameter 
space. 
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,„ nil 



FIG. 5: In this group of figures, we look at the shapes near the double squeezed limit: we consider the 
case where ks = k^ = k\i and the tetrahedron is a planar quadrangle. We plot T; OCj eq , Ti oc \ and T/ OC 2, 
respectively, as functions of k^jk\ and ku/ki. Notice that the sign of T} oc , eq varies non trivially over the 
domain. In the squeezed limit, at (k^jk\ — l,/c 1 4/fc 1 = 1) where /c 2 — > 0, and in the double-squeezed limit, 
&3 = ki — > 0, Ti OCi eg blows up while the local shapes are finite (see the main text for additional comments). 
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the tetrahedron reduces to a planar quadrangle, its sign non-trivially varies over the parameter 
space, in contrast with all shapes studied in [78] for instance. This is all the more interesting 
as this planar limit corresponds to the situation effectively probed by small (angular) scale CMB 



experiments. This gives hope to be able to test the consistency relation (94). 



V. CONCLUSIONS 



In this paper, we have analyzed a scenario to convert entropic into curvature perturbations in 
the context of multifield brane inflation, and multifield DBI inflation in particular. Considering a 
standard single-field driven inflationary phase ending by tachyonic instability, we have shown how 
the perturbations of light fields, parametrizing the angular position of the brane in the throat, can 
generate curvature perturbation by spatially modulating the duration of inflation. The entropic 
transfer generated by this mechanism is more efficient in the DBI than in the slow-roll regime, 
being enhanced by the inverse of the sound speed at horizon crossing and at the end of inflation. 

We have investigated the non-Gaussianities created in such models. Single-field DBI inflation 
is known to produce a substantial amount of non-Gaussianities of equilateral type and previous 
works on multifield DBI inflation showed that the entropic perturbations reduce the amplitude 
fS L of these non-Gaussianities. However, as soon as light scalar fields other than the inflaton are 
present, non-Gaussianities of another shape, namely local non-Gaussianities, can be present. In 
our model, they arise because of the nonlinear relation between the curvature perturbation and 
the angular ones. This local contribution to the non-linearity parameter /jvl can be large and 
even saturate the observational bounds for some specific scenarios. 

We have also calculated tnl and gNL, which are trispectrum parameters similar to for the 
bispectrum, in that they describe the nonlinearities generated classically outside the horizon. We 
have shown that qnl is related to tnl by a simple parameter describing the angular position of 
the mobile brane, implying g^i < j| tnl • When the curvature perturbation is mostly of entropic 
origin, t^l = (f/jvi,) 2 ' easu y reaching the expected Planck sensitivity \tnl\ ~ 560 [75] as soon as 
JNL ~ 25. The trispectrum also contains a part that is directly related to the trispectra of the field 
perturbations. One can in principle extract the entropic transfer function from its complicated 
momentum dependence, as well as determine the sound speed when the observable modes cross 
the horizon by combining it with measurements of f^ L - 

Finally, due to the presence of light scalar fields, other than the inflaton, and with non standard 
kinetic terms, the primordial trispectrum acquires a particular momentum dependent component 
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whose amplitude (94) is given by the product f^L- Moreover, this consistency relation is valid 
in every multifield DBI model, whatever the brane's trajectory or the mechanism to convert entropic 
perturbations into the curvature perturbation. It thus represents an interesting observational 
signature of multifield DBI inflation. We have represented the corresponding form factor in different 
limits and it turned out to display important differences with the shapes associated to ttvl an d gNL, 
such as its sign varying over the tetrahedron's parameter space. Overall, should future experiments 
come to detect both local and equilateral non-Gaussianities in the primordial bispectrum, the 
trispectrum may well help to to confirm or exclude this kind of scenarios. 

Note added: On the day this work was submitted, the paper [B6 appeared in the arXiv, which 
completes the calculation of the trispectrum coming from the four-point functions of the field 
perturbations in multifield DBI inflation. 
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Appendix 1 - 5N formulae beyond leading order 



Using the expressions of (29) and ( e (30) together with the definitions of rj (12) and s 



(13), one obtains the coefficients of the SN expansion without restricting to leading order in these 



slowly- varying parameters and their time derivatives: 



N„ 



1 



1 



d>L l 



(98) 
(99) 



Nrr 



N.. 



7] + S 1 



4ec, M p 2 
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^2e~c~ s M pX 4ec s M p 2 L 



(100) 
(101) 
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N ™° ~ 2 5/2( eCs )3/2 M 3 [ H + H + ^ 

n sss =-(^£ +3 fa+')^ w 

+ 2 5/2( eCs) 3/2 M 3 ^ + # - W + S ) 



(102) 

) 

v/2l^M p i ' " 4ec s M*~ 

(103) 



The derivatives of (f) e (E*) (|26j) are evaluated at the background value S = b*c s *0 and their 
explicit expressions read 



/3 

tan(a)— , (104) 



cos d (a) \Cs* / 0c 
(3) _ 3tan(a) / /? \ 3 1 



cos 4 (a) \ c s*/ 0c 
Appendix 2 - Higher order loop contributions 



(105) 
(106) 



Given we encountered large entropic derivatives in the 5N formalism, one must ensure that we 
are in a simple perturbative regime in which one can safely neglect higher order loop corrections 
(see |87[ [88] for discussions on the viability of a perturbative expansion and (89 J for a diagrammatic 
approach to loop effects). The integral over the loop momenta give rise to logarithmic infrared 
divergences ln(fcL) where L is the large scale cut off. In the following we assume that one can take 



ln(fc-L) ~ 1 as done in \90\ l9T| [92] . Concentrating on the 8N expansion (24) to second order for 
simplicity, and neglecting for the moment the non-Gaussianities of the fields at horizon crossing, 
the dominant loop corrections to the power spectrum, bispectrum and trispectrum are |91[ [92] 

flloop 



Pj " N AB N AB , 

ptree " (iVc jVc) 2 



rl loop 

JNL ~ ft fAT-AT-\3 'O I 108 ) 



pi loop 

5N AB N BC N A c, 
6 (N D N D y 
Hoop N ab N bc N C dN ad 



Below we discuss the two limiting cases of a large and small entropic transfer. 
• Large entropic transfer. 
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For a large entropic transfer, one obtains 



p ll00p / AT \ 2 / a \ 2 

fC = f ^££ \ -P = ( - 

ptree \m ) C I 5 



C 



.1 loop _ 5 / 7V SS \ _ 36 / > oc \ 3 

./A/7, — o A79 ' C — or Uj 



n = ( ) r (110) 

3 



rlir = (^) lp < = (^Ai) 4 n. (112) 

Therefore, given that Vt ~ 10" 10 and that observations already show that the level of non- 
Gaussianities is relatively small, = 0(100), this demonstrates that loop corrections can 
be neglected, and clearly the argument is valid for every scenario in which one Gaussian field 
only is responsible for the curvature perturbation |93j . 

In general, the field non-Gaussianities lead to additional loop contributions, for instance 
to Inl- This was shown to be completely negligible for slow-roll models [94] but in DBI 
inflation where intrinsic non-Gaussianities are large, one has to determine if this remains 



true. In particular, there is a contribution to the primordial trispectrum (46) from terms 
of the form 



N a NbcN de (Q a (Mi)(Q B *Q c )(k 2 )(Q D *Q*)(k 3 )> , (113) 



where the symbol * denotes a convolution product. The most dangerous terms in (113) 
involve the entropic derivatives and, as to leading order there is no purely entropic three 
point function, one is led to consider 

Af CT ^ 2 s (Q (7 (k 1 )(Q^Q s )(k 2 )(Q s *g s )(k3)). (114) 



From the definition of /jvx (48) as well as the result (37) for the three point functions of the 



fields, its contribution is of order 

(iV^-^g-^)/( ff ,X 4 (l + ^) 2 ). ("5) 
For a large entropic transfer, one therefore obtains 

fNL° P ^ fN q L(f l NL) 2 r C , (116) 



where we have used (52) and (53). Again, the observational bounds = 0(100) as 

well as the normalization V( ~ 10" 10 show that such a loop-correction is negligible and one 
can verify that the same conclusion applies to other higher order corrections. 
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Small entropic transfer 



When the entropic transfer is inefficient - T^ s <1- the power spectrum ( 42 ) as well as the 



equilateral non-Gaussianities (52) remain the same as in the single-field case, so that the 
multifield aspects we are taking into account can not cure the model if it is under pressure 



at these levels. As for the local non-Gaussianity parameter, from (53) one obtains 

5. 



lot: 
NL 



_ T 2 Nbs 

6 as m ' 



TL < l 



(117) 



Its suppression by the small entropic transfer gives few hope to have a significant 
rendering the situation effectively indistinguishable from the single-field case for these ob- 
servables, unless we have an extremely large ratio N ss /N%. For instance, T^ s ~ 10 -2 and 
N ss /N% ~ 10 3 give ffil ~ 10. However, if N ss /N% is made this large, Eqs. ( jlOg) and ( fl09| ) 
give 



A loop 
J NL 



1 loop 
T NL 



5 fN, 



6 \N$ J 

m 



(118) 
(119) 



and loop effects become important, which requires a more refined study. 



[1] E. Komatsu et al. (2009), 0902.4759. 

[2] P. Creminelli and M. Zaldarriaga, JCAP 0410, 006 (2004), astro-ph/0407059. 

[3] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 180, 330 (2009), 0803.0547. 

[4] D. Baumann, A. Dymarsky, I. R. Klcbanov, and L. McAllister, JCAP 0801, 024 (2008), 0706.0360. 

[5] D. Baumann, A. Dymarsky, I. R. Klebanov, L. McAllister, and P. J. Steinhardt, Phys. Rev. Lett. 99, 
141601 (2007), 0705.3837. 

[6] E. Silverstein and D. Tong, Phys. Rev. D70, 103505 (2004), hep-th/0310221. 

[7] M. Alishahiha, E. Silverstein, and D. Tong, Phys. Rev. D70, 123505 (2004), hep-th/0404084. 

[8] R. Bean, S. E. Shandera, S. H. Henry Tye, and J. Xu, JCAP 0705, 004 (2007), hep-th/0702107. 

[9] H. V. Peiris, D. Baumann, B. Friedman, and A. Cooray, Phys. Rev. D76, 103517 (2007), 0706.1240. 

[10] A. Kehagias and E. Kiritsis, JHEP 11, 022 (1999), hep-th/9910174. 

[11] D. A. Steer and M. F. Parry, Int. J. Thcor. Phys. 41, 2255 (2002), hep-th/0201121. 

[12] P. Brax and D. A. Steer, Phys. Rev. D66, 061501 (2002). 

[13] D. Easson, R. Gregory, G. Tasinato, and I. Zavala, JHEP 04, 026 (2007), hep-th/0701252. 

[14] C. Gordon, D. Wands, B. A. Bassett, and R. Maartens, Phys. Rev. D63, 023506 (2001), astro- 
ph/0009131. 



35 



[15] D. Langlois, S. Renaux-Petel, D. A. Steer, and T. Tanaka, Phys. Rev. Lett. 101, 061301 (2008), 
0804.3139. 

[16] D. Langlois, S. Renaux-Petel, D. A. Steer, and T. Tanaka, Phys. Rev. D78, 063523 (2008), 0806.0336. 

[17] F. Arroja, S. Mizuno, and K. Koyama, JCAP 0808, 015 (2008), 0806.0619. 

[18] D. Langlois, S. Renaux-Petel, and D. A. Steer, JCAP 0904, 021 (2009), 0902.2941. 

[19] X. Gao and B. Hu (2009), 0903.1920. 

[20] X. Gao and F. Xu (2009), 0905.0405. 

[21] S. Mizuno, F. Arroja, K. Koyama, and T. Tanaka (2009), 0905.4557. 

[22] X. Gao, M. Li, and C. Lin (2009), 0906.1345. 

[23] D. H. Lyth and A. Riotto, Phys. Rev. Lett. 97, 121301 (2006), astro-ph/0607326. 

[24] A. Sen, JHEP 07, 065 (2002), hep-th/0203265. 

[25] H.-Y. Chen, J.-O. Gong, and G. Shiu, JHEP 09, 011 (2008), 0807.1927. 

[26] L. Lcblond and S. Shandcra, JCAP 0701, 009 (2007), hcp-th/06 10321. 

[27] D. Baumann and L. McAllister, Phys. Rev. D75, 123508 (2007), hcp-th/0610285. 

[28] J. E. Lidsey and I. Huston, JCAP 0707, 002 (2007), 0705.0240. 

[29] S. B. Giddings, S. Kachru, and J. Polchinski, Phys. Rev. D66, 106006 (2002), hcp-th/0105097. 

[30] I. R. Klcbanov and M. J. Strasslcr, JHEP 08, 052 (2000), hcp-th/0007191. 

[31] D. Baumann and L. McAllister (2009), 0901.0265. 

[32] D. Langlois and S. Renaux-Petel, JCAP 0804, 017 (2008), 0801.1085. 

[33] V. F. Mukhanov, H. A. Fcldman, and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 

[34] D. Langlois (2004), hcp-th/0405053. 

[35] O. DeWolfe, S. Kachru, and H. L. Vcrlinde, JHEP 05, 017 (2004), hcp-th/0403123. 

[36] O. DeWolfe, L. McAllister, G. Shiu, and B. Underwood, JHEP 09, 121 (2007), hcp-th/0703088. 

[37] E. Pajer, JCAP 0804, 031 (2008), 0802.2916. 

[38] F. Bernardeau and J.-P. Uzan, Phys. Rev. D66, 103506 (2002), hcp-ph/0207295. 

[39] M. Sasaki, Prog. Thcor. Phys. 120, 159 (2008), 0805.0974. 

[40] A. Naruko and M. Sasaki, Prog. Thcor. Phys. 121, 193 (2009), 0807.0180. 

[41] C. T. Byrnes, K.-Y. Choi, and L. M. H. Hall, JCAP 0902, 017 (2009), 0812.0807. 

[42] Q.-G. Huang, JCAP 0905, 005 (2009), 0903.1542. 

[43] A. A. Starobinsky, JETP Lett. 42, 152 (1985). 

[44] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996), astro-ph/9507001. 

[45] M. Sasaki and T. Tanaka, Prog. Thcor. Phys. 99, 763 (1998), gr-qc/9801017. 

[46] D. H. Lyth, K. A. Malik, and M. Sasaki, JCAP 0505, 004 (2005), astro-ph/0411220. 

[47] D. H. Lyth and Y. Rodriguez, Phys. Rev. Lett. 95, 121302 (2005), astro-ph/0504045. 

[48] D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddlc, Phys. Rev. D62, 043527 (2000), astro- 
ph/0003278. 

[49] D. H. Lyth, JCAP 0511, 006 (2005), astro-ph/0510443. 



36 



[50] N. Barnaby and J. M. Clinc, Phys. Rev. D75, 086004 (2007), astro-ph/ 06 11750. 

[51] K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki, and A. Vaihkonen, JHEP 08, 084 (2005), 
hep-th/0502185. 

[52] C. T. Byrnes, M. Sasaki, and D. Wands, Phys. Rev. D74, 123519 (2006), astro-ph/0611075. 

[53] J. M. Maldacena, JHEP 05, 013 (2003), astro-ph/0210603. 

[54] D. Seery and J. E. Lidscy, JCAP 0509, Oil (2005), astroph/0506056. 

[55] X. Chen, M.-x. Huang, S. Kachru, and G. Shiu, JCAP 0701, 002 (2007), hcp-th/0605045. 

[56] X. Chen, Phys. Rev. D72, 123518 (2005), astro-ph/0507053. 

[57] D. Wands, N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D66, 043520 (2002), astro-ph/0205253. 

[58] S. Kecskcmcti, J. Maiden, G. Shiu, and B. Underwood, JHEP 09, 076 (2006), hcp-th/0605189. 

[59] L. E. Allen, S. Gupta, and D. Wands, JCAP 0601, 006 (2006), astro-ph/0509719. 

[60] Q.-G. Huang (2009), 0904.2649. 

[61] F. Bernardeau and J.-P. Uzan, Phys. Rev. D67, 121301 (2003), astro-ph/0209330. 

[62] G. Dvali, A. Gruzinov, and M. Zaldarriaga, Phys. Rev. D69, 023505 (2004), astro-ph/0303591. 

[63] L. Kofman (2003), astro-ph/0303614. 

[64] M. Zaldarriaga, Phys. Rev. D69, 043508 (2004), astro-ph/0306006. 

[65] F. Bernardeau, L. Kofman, and J.-P. Uzan, Phys. Rev. D70, 083004 (2004), astro-ph/0403315. 

[66] M. P. Salem, Phys. Rev. D72, 123516 (2005), astro-ph/0511146. 

[67] L. Alabidi and D. Lyth, JCAP 0608, 006 (2006), astro-ph/0604569. 

[68] F. Bernardeau and T. Brunicr, Phys. Rev. D76, 043526 (2007), 0705.2501. 

[69] B. Dutta, L. Lcblond, and J. Kumar, Phys. Rev. D78, 083522 (2008), 0805.1229. 

[70] D. Babich, P. Crcminclli, and M. Zaldarriaga, JCAP 0408, 009 (2004), astro-ph/0405356. 

[71] M. LoVerde, A. Miller, S. Shandcra, and L. Verde, JCAP 0804, 014 (2008), 0711.4126. 

[72] E. Sefusatti, M. Liguori, A. P. S. Yadav, M. G. Jackson, and E. Pajer (2009), 0906.0232. 

[73] S. Bird, H. V. Pciris, and D. Baumann (2009), 0905.2412. 

[74] A. R. Frey, A. Mazumdar, and R. C. Myers, Phys. Rev. D73, 026003 (2006), hcp-th/0508139. 

[75] N. Kogo and E. Komatsu, Phys. Rev. D73, 083007 (2006), astro-ph/0602099. 

[76] A. Cooray, C. Li, and A. Melchiorri, Phys. Rev. D77, 103506 (2008), 0801.3463. 

[77] D. Jcong and E. Komatsu (2009), 0904.0497. 

[78] X. Chen, B. Hu, M.-x. Huang, G. Shiu, and Y. Wang (2009), 0905.3494. 

[79] T. Okamoto and W. Hu, Phys. Rev. D66, 063008 (2002), astro-ph/0206155. 

[80] M. Sasaki, J. Valiviita, and D. Wands, Phys. Rev. D74, 103003 (2006), astro-ph/0607627. 

[81] J.-L. Lchners and S. Rcnaux-Pctel (2009), 0906.0530. 

[82] C. T. Byrnes and G. Tasinato (2009), 0906.0767. 

[83] D. Seery, M. S. Sloth, and F. Vcrnizzi, JCAP 0903, 018 (2009), 0811.3934. 

[84] F. Arroja, S. Mizuno, K. Koyama, and T. Tanaka (2009), 0905.3641. 

[85] X. Chen, M.-x. Huang, and G. Shiu, Phys. Rev. D74, 121301 (2006), hcp-th/0610235. 



37 



[86] S. Mizuno, F. Arroja, and K. Koyama (2009), 0907.2439. 

[87] C. Cheung, P. Creminelli, A. L. Fitzpatrick, J. Kaplan, and L. Senatore, JHEP 03, 014 (2008), 
0709.0293. 

[88] L. Lcblond and S. Shandera, JCAP 0808, 007 (2008), 0802.2290. 

[89] C. T. Byrnes, K. Koyama, M. Sasaki, and D. Wands, JCAP 0711, 027 (2007), 0705.4096. 

[90] L. Boubekeur and D. H. Lyth, Phys. Rev. D73, 021301 (2006), astroph/0504046. 

[91] H. R. S. Cogollo, Y. Rodriguez, and C. A. Valenzuela-Toledo, JCAP 0808, 029 (2008), 0806.1546. 

[92] Y. Rodriguez and C. A. Valenzuela-Toledo (2008), 0811.4092. 

[93] D. H. Lyth, JCAP 0712, 016 (2007), 0707.0361. 

[94] I. Zaballa, Y. Rodriguez, and D. H. Lyth, JCAP 0606, 013 (2006), astro-ph/0603534. 
[95] X. Chen, Phys. Rev. D71, 063506 (2005), hep-th/0408084. 
[96] X. Chen, JHEP 08, 045 (2005), hep-th/0501184. 

[97] R. Bean, X. Chen, H. Peiris, and J. Xu, Phys. Rev. D77, 023527 (2008), 0710.1812. 

[98] X. Chen, S. Sarangi, S. H. Henry Tye, and J. Xu, JCAP 0611, 015 (2006), hep-th/0608082. 

[99] J. Khoury and F. Piazza (2008), 0811.3633. 
[100] S. Renaux-Petel and G. Tasinato, JCAP 0901, 012 (2009), 0810.2405. 
[101] D. Langlois and F. Vernizzi, JCAP 0702, 017 (2007), astro-ph/0610064. 
[102] Y.-F. Cai and W. Xue (2008), 0809.4134. 

[103] Y.-F. Cai and H.-Y. Xia, Phys. Lett. B677, 226 (2009), 0904.0062. 
[104] L. Senatore, K. M. Smith, and M. Zaldarriaga (2009), 0905.3746. 

[105] We are considering the so-called ultra-violet model of DBI inflation, in which the mobile D3-brane falls 
towards the tip of the throat, but consistent infrared models [95j [W] can be constructed, as explained 
in [37]. 

[106] This bound was first discussed in [98] in the context of eternal D-brane inflation. 

[107] We use the symbol p to evade any ambiguities with the commonly used variable r in the Klebanov- 
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[108] See [55] for models where the speed of sound is rapidly varying. 

[109] Notice that for the KS throat, its metric is often described with another radial variable, r [3D], such 
that p ~ \/g s Mh t ipl s T near the tip, where M ^> 1 is a flux integer. With this variable, the tachyon 
appears - forgetting about the angular direction - for r ~ l/\/g s M, where the warped string length 
htiph does not appear. This field however is not canonically normalized and, had we used this variable, 
the warped string length hti p l s would reappear in the amplitude of the corresponding perturbation. 
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[Ill] Note that strictly speaking there exists also a purely adiabatic contribution to /Jyf,, proportional to 



era- ( 100 1 . This equals ^(i+j^* )2 , so it is suppressed by both slow- varying parameters and by the 
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cntropic transfer and is hence unobservably small. Similar, purely adiabatic, contributions are present 

in the non-linearity parameters of the trispectrum and we will neglect them as well. 

[112] Note that this could be realized in the different context of multiple branes inflation [1021 1103] , 

[113] We thank Eiichiro Komatsu for pointing this to us. Note also the recent paper |104j which demonstrates 

that the analysis of the shapes of non-Gaussianities requires careful handling. 

[114] A more severe bound can even be derived if one requires that the four-dimensional energy-density p e 

be less than the local string energy density {hti p m s ) A . With the Friedmann equation 3M^H^ = p e , this 

(M \ 1/2 

gives indeed h t i p m s /H e > I -jf 1 ) > 140 where the last inequality follows from the non detection of 
tensor modes. As this is only an order of magnitude bound, we do not consider that it rules out the 
model. 



